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HIGHLIGHTS 


►  Flooding  characteristics  for  cathode  GDL  of  direct  humidified  PEMFC. 

►  Co-flow  and  counter-flow  of  fuel  and  oxidant  for  pressure  conditions  of  1  and  3  atm. 

►  Numerical  analysis  results  showed  that  flooding  area  was  limited. 


ARTICLE  INFO 


ABSTRACT 


Article  history: 

Received  14  October  2011 
Received  in  revised  form 
20  December  2011 
Accepted  26  April  2012 
Available  online  5  May  2012 


Keywords: 

Proton  exchange  membrane  fuel  cell 

Direct  humidified  cell 

Flooding 

Numerical  analysis 


In  proton  exchange  membrane  fuel  cell  (PEMFC)  stack,  the  water  management  in  the  cell  is  required  in 
order  to  ensure  high  cell  performance.  The  water  dry-out  interferes  with  the  proton  migration  in  the 
electrode  and  electrolyte,  and  the  water  flooding  interferes  with  the  reactant  gas  transfer  in  the  gas 
diffusion  layer,  consequently  these  phenomena  affect  the  cell  performance.  One  of  the  approaches  to 
prevent  the  dry-out,  a  direct  humidified  cell  with  the  channel  to  supply  liquid  phase  water  has  been 
considered.  In  this  cell  configuration,  sufficient  water  will  be  supplied  to  the  anode,  however  the 
migrated  water  to  the  cathode  from  the  anode  shall  be  removed  immediately  in  order  to  prevent  the 
flooding.  In  this  paper,  numerical  analysis  that  includes  the  behavior  of  the  electrode  reaction,  heat 
transfer,  gas  flow,  and  evaporation  and  condensation  of  water  was  conducted  regarding  a  direct 
humidified  model  cell.  The  co-flow  and  counter-flow  of  fuel  and  oxidant  gases  were  considered  for  the 
pressure  conditions  of  1  atm  and  3  atm.  The  water  distribution  in  the  vapor  phase  and  liquid  phase  at  the 
oxidant  gas  and  cathode  were  calculated,  and  as  a  result  of  those,  the  flooding  characteristics  for  each 
condition  were  evaluated. 

©  2012  Published  by  Elsevier  B.V. 


1.  Introduction 

In  proton  exchange  membrane  fuel  cell  (PEMFC),  it  is  well 
known  that  appropriate  humidity  in  the  cell  is  required  in  order  to 
ensure  high  cell  performance.  The  dry-out  caused  by  low  humidity 
interferes  with  the  proton  migration  in  the  membrane  electrolyte 
assembly  (MEA),  and  the  flooding  caused  by  excessive  humidity 
interferes  with  the  reactant  gas  transfer  in  the  gas  diffusion  layer 
(GDL).  Therefore,  studies  regarding  the  liquid  and  vapor  phase 
water  management  for  the  MEA,  GDL,  and  gas  channel  have  been 
conducted  in  these  years. 

For  the  purpose  of  evaluating  the  cathode  channel  flooding, 
numerous  studies  have  been  conducted.  Shimpulee  et  al.  [1]  have 
conducted  numerical  study  of  the  flooding  in  the  GDL  that  includes 
a  micro  layer  porous  electrode  using  computational  fluid  dynamics 
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(CFD)  tool,  and  discussed  about  the  cell  performance  reduction 
caused  by  various  filling  level  of  liquid  water  for  the  electrode  pore. 
Le  et  al.  [2]  have  conducted  numerical  study  of  the  flooding  in  the 
gas  channel  and  GDL  using  CFD  tool,  which  includes  the  liquid 
phase  water  migration,  and  investigated  the  liquid  water  existence 
in  the  serpentine  gas  channels.  Berg  et  al.  [3]  have  conducted 
numerical  study  of  the  liquid  and  vapor  phase  water  transfer 
characteristics  in  the  cathode  channel  using  one  dimensional 
model  of  the  co-flow  and  counter-flow  for  the  fuel  and  oxidant.  Cho 
et  al.  [4]  have  conducted  coupling  study  of  the  experiment  and  CFD 
calculation,  and  estimated  the  flooding  region  with  calculation 
based  on  the  cell  surface  temperature  that  have  been  measured  by 
experiment.  Chen  et  al.  [5]  have  conducted  dynamic  simulation 
study  of  the  cell  and  humidifier,  and  calculated  dynamic  response 
of  the  vapor  transfer  in  the  membrane. 

On  the  other  hand,  for  the  purpose  of  evaluating  dry-out,  Lister 
et  al.  [6]  have  conducted  experimental  study,  and  dry-out 
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characteristics  have  been  observed  by  measuring  cell  resistance 
with  dry  gas  flow. 

A  direct  humidified  cell  in  which  the  liquid  phase  water  is 
supplied  has  been  considered  for  an  approach  to  prevent  the  dry¬ 
out.  In  this  cell  configuration,  the  water  channel  that  was  con¬ 
tacted  to  the  fuel  channel  will  supply  sufficient  water  to  the  anode 
GDL,  and  the  dry-out  will  be  prevented  effectively.  However,  the 
migrated  water  to  the  cathode  GDL  and  oxidant  channel  from  the 
anode  GDL  shall  be  removed  immediately  by  the  evaporation  with 
oxidant  gas  flow  in  order  to  prevent  the  flooding. 

In  this  paper,  numerical  analysis  that  includes  the  behavior  of  the 
electrode  reaction,  heat  transfer,  gas  flow,  and  evaporation  and 
condensation  of  water  was  conducted  regarding  a  direct  humidified 
model  cell.  The  co-flow  and  counter-flow  of  fuel  and  oxidant  gases 
were  considered  for  the  pressure  conditions  of  1  atm  and  3  atm.  The 
flooding  areas  regarding  each  condition  were  estimated,  and  as  a  result 
of  those  the  effect  of  flooding  on  the  cell  performance  was  evaluated. 

2.  Analysis  method 

The  analysis  features  about  the  heat  transfer  and  gas  flow  in  the 
stack  are  the  same  as  those  of  the  code  for  MCFC  [8],  and  the 
functions  of  the  water  evaporation  and  condensation  have  been 
added  based  on  the  code  for  PAFC  [7],  and  function  of  the  electrode 
reaction  has  been  modified  to  simulate  the  PEMFC. 

Regarding  the  electrode  reaction,  the  following  V—i  character¬ 
istic  equation  [9]  was  adopted.  This  equation  has  been  derived  from 
the  experiments  of  the  single  cell  using  50.56  cm2  “Nation” 
membrane  for  the  conditions  of  current  density  54-323  mA  cm-2, 
temperature  323-358  K,  hydrogen  partial  pressure  1.90-3.06  atm, 
and  oxygen  partial  pressure  0.64-3.06  atm. 

V  =  1.229-0.85  X  10  3  x  (T  -  298.15)  +  4.3085  x  10“5  xT 
x  [ln(xH2  xx°f)+ln(p15)]  -0.9514  +  3.12  x  103  x  T 
-  1.87  x  10~4  x  T  x  ln(s  x  i)  +  7.4  x  10~5  xTx  ln(c0J 
-s  xi  x  (0.01605  -  3.5  x  10-5  x  T  +  8.0  x  10-5  x  s  xi) 

0) 

Here,  V  is  cell  voltage  V,  T  is  temperature  K,  x  is  molar  fraction,  P 
is  pressure  atm,  c  is  molar  concentration  mol  cm-2,  s  is  cell  area 
50.56  cm2,  and  i  is  current  density  A  cm-2. 

The  water  evaporation  and  condensation  rates  were  assumed  to 
be  infinite  rates  based  on  the  saturated  condition  for  the  liquid 
phase  water  at  the  GDL  and  the  vapor  phase  water  at  the  oxidant  gas. 

The  phosphoric  acid  evaporation  and  condensation  rates  in  the 
PAFC  had  been  calculated  based  on  the  activity  energy  change  of 
the  two  components  system  consists  of  the  water  and  phosphorous 
pentoxide,  and  the  evaporation  and  condensation  rates  were 
derived  with  sufficient  accuracy  [7].  Where,  the  activity  energy 
change  was  derived  based  on  the  vapor  liquid  equilibrium  condi¬ 
tion  of  the  water  and  phosphorous  pentoxide.  Therefore,  the  water 
evaporation  and  condensation  rates  in  this  paper,  in  which  the 
equilibrium  condition  were  assumed,  were  also  considered  to  have 
sufficient  accuracy. 

In  addition  to  the  above  assumption,  the  following  conditions 
were  assumed. 

(a)  The  anode  GDL  was  assumed  to  be  filled  with  sufficient  water 
from  the  water  channel. 

(b)  The  amount  of  transferred  water  from  anode  to  cathode  by  the 
drag  of  proton  and  diffusion  was  assumed  to  be  1.4  times  the 
amount  of  composed  by  the  electrode  reaction. 
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Fig.  1.  Model  cell  configuration. 


(c)  The  composed  water  at  the  cathode  was  assumed  to  be  the 
liquid  phase. 

(d)  The  amount  of  evaporating  water  from  the  cathode  GDL  to 
oxidant  gas  was  assumed  to  be  less  than  or  equal  to  the  amount 
of  water  that  was  migrated  from  the  anode  GDL. 

The  assumption  (b)  was  based  on  the  results  of  the  amount  of 
migrated  water  from  the  anode  to  cathode  by  the  drag  of  proton 
and  diffusion,  which  was  almost  1.4  times  the  composed  water 
with  the  electrode  reaction  [10].  Where,  it  was  derived  regarding 
various  current  density  conditions  from  0.05  to  1.7  A  cm-2  for 
conditions  of  90  °C  and  2  atm. 

3.  Analysis  model  and  conditions 

The  outline  of  model  cell  is  shown  in  Fig.  1,  and  specifications  are 
shown  in  Table  1.  The  model  cell  consists  of  the  anode  GDL,  MEA, 
cathode  GDL,  anode  separator  with  fuel  gas  channel,  cathode 
separator  with  oxidant  gas  channel,  and  the  separator  including 


Table  1 

Specifications  for  the  model  cell. 


Cell  size 

170  x  170  mm 

Separator  thickness 

1.5  mm 

Channel  width  and  height 

1.0  x  0.5  mm 

GDL  thickness 

0.3  mm 

Separator  thermal  conductivity 

139.6  W  m”1  K1  (in  plane  direction) 
4.7  W  m-1  K  1  (thickness  direction) 

H.  Hirata  et  al.  /  Journal  of  Power  Sources  220  (2012)  383-390 


385 


water  channel.  The  each  separator  has  18  serpentine  channels  that 
turn  2  times.  In  the  case  of  co-flow  configuration,  the  fuel,  water, 
and  oxidant  were  supplied  at  the  left  bottom  corner  and  discharged 
at  the  right  top  corner  in  Fig.  1.  In  the  case  of  counter-flow 
configuration,  the  fuel  and  water  were  supplied  at  the  left  bottom 
corner  and  discharged  at  the  right  top  corner,  and  the  oxidant  was 
supplied  at  the  right  top  corner  and  discharged  at  the  left  bottom 
corner  in  Fig.  1.  The  repeating  thermal  boundary  condition  that  the 
water  and  oxidant  separator  surfaces  had  the  same  temperature 
and  temperature  gradient  was  adapted,  as  a  result  of  this;  the  cell 
model  represented  a  stack  consists  of  infinity  stacking  cells.  The 
adiabatic  condition  was  adopted  in  each  side  of  the  cell.  Analysis 
was  conducted  regarding  the  co-flow  and  counter-flow  configura¬ 
tions  for  the  pressure  conditions  of  1  and  3  atm. 

The  other  conditions  are  shown  in  Table  2.  Flere,  the  current 
density  condition  of  350  mA  cm'2  was  selected  to  correspond  with 
the  equation  (1 ).  As  the  main  purpose  of  this  paper  was  to  evaluate 
the  effect  of  the  water  evaporation  and  condensation  on  the 
flooding  area,  a  comparative  study  for  such  relatively  low  current 
density  condition  was  attempted. 

4.  Analysis  results 

4.1  Cell  voltage  and  cell  average  temperature 

The  cell  voltage  and  cell  average  temperature  for  each  analysis 
case  are  shown  in  Table  3.  The  cell  voltage  was  increased  about  7.7% 
and  the  cell  average  temperature  decreased  1.8  K,  according  to  the 
pressure  increase  from  1  atm  to  3  atm  for  each  gas  flow  configu¬ 
ration.  The  effect  of  the  gas  flow  configurations  on  the  cell  voltage 
and  cell  average  temperature  was  few. 


4.2.  Current  density,  cell  temperature,  water  evaporation  rate,  and 
water  evaporation  margin 


The  current  density  distribution,  cell  temperature  distribution, 
water  evaporation  rate  distribution  on  cathode  GDL,  and  water 
evaporation  margin  distribution  on  cathode  GDL  for  the  conditions 
of  1  atm  and  co-flow  are  shown  in  Figs.  2-5,  respectively. 

Flere,  net  water  evaporation  or  condensation  rate  was  derived  as 
infinite  value  based  on  the  water  equilibrium  conditions  as 
mentioned  in  Section  2.  Flowever,  as  a  result  of  combination  with 
the  oxidant  gas  flow  rate,  the  following  apparent  water  evaporation 
rate  was  derived. 


t  _  P h20  v  Qo 
•*a  P  A 


(2) 


Flere,  pH2o  is  saturated  water  vapor  pressure  atm,  P  is  total 
pressure  atm,  Qo  is  oxidant  molar  flow  rate  mol  s'1,  A  is  cathode 
GDL  area  cm2,  and  Ja  is  apparent  water  evaporation  rate 
mol  s'1  cm'2.  As  shown  in  equation  (2),  the  apparent  water 
evaporation  rate  increased  according  to  the  oxidant  gas  flow  rate 
increase,  and  decreased  according  to  the  total  pressure  increase. 

The  water  evaporation  margin  on  cathode  GDL  was  defined  as 
follows; 


Table  2 

Analysis  conditions. 


Pressure 

1  atm,  3  atm 

Gas  inlet  temperature 

80  °C 

Average  current  density 

350  mA  cm-2 

Fuel  utilization  (Dry  H2) 

70% 

Oxidant  utilization  (Dry  air) 

40% 

Gas  flow  configuration 

Co-flow,  Counter-flow 

Table  3 

Analysis  results. 


Analysis  conditions  Analysis  results 


Pressure  (atm)  Gas  flow  Cell  voltage  (V)  Average  cell  temperature  (K) 


1.0 

Co-flow 

0.6813 

358.1 

1.0 

Counter-flow 

0.6814 

358.2 

3.0 

Co-flow 

0.7336 

356.4 

3.0 

Counter-flow 

0.7336 

356.4 

Gn  —  Q— >o  —  Cq— >c  (3) 

Here,  Cm  is  water  evaporation  margin  mol  s'1  cm'2,  Cc^0  is  the 
amount  of  transferred  water  from  the  cathode  GDL  to  oxidant  gas 
mol  s'1  cm'2,  and  Cc^0  is  the  amount  of  migrated  water  from  the 
anode  GDL  to  cathode  GDL  mol  s'1  cm'2.  In  the  case  of  Cm  is  plus, 
the  amount  of  migrated  water  from  the  anode  to  cathode  will  be 
fully  transferred  to  the  oxidant  gas  and  the  possibility  of  flooding 
will  be  low.  On  the  other  hand,  in  the  case  of  Cm  is  minus,  the 
amount  of  migrated  water  from  the  anode  to  cathode  will  not  be 
fully  transferred  to  the  oxidant  gas  and  the  possibility  of  flooding 
will  be  high. 

The  current  density  was  the  largest  at  the  inlet  of  oxidant  and 
fuel,  and  it  decreased  16%  through  the  serpentine  channel  for  the 
outlet  of  oxidant  and  fuel  in  Fig.  2.  Here,  the  current  density 
distribution  was  determined  by  equation  (1)  according  to  the 
pressure  condition  and  the  temperature  distribution  and  molar 
fraction  distributions  of  H2  and  02.  The  cell  temperature  was  the 
lowest  at  the  inlet  of  oxidant  and  fuel,  and  it  increased  0.84% 
through  the  serpentine  channel  for  the  outlet  of  oxidant  and  fuel  in 
Fig.  3,  because  the  oxidant  gas  transferred  the  increased  enthalpy 
by  the  electrode  reaction  from  the  oxidant  inlet  to  the  oxidant 
outlet. 

The  apparent  water  evaporation  rate  was  the  largest  at  the  inlet 
of  oxidant  and  fuel,  and  it  decreased  18%  through  the  serpentine 
channel  for  the  outlet  of  oxidant  and  fuel  in  Fig.  4.  The  water  partial 


Fuel 


Fuel 

Fig.  2.  Current  density  (A  cm  2),  1  atm,  co-flow. 
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Fig.  3.  Cell  temperature  (K),  1  atm,  co-flow. 


Fig.  5.  H20  evaporation  margin  (mol  s  1  cm  2),  1  atm,  co-flow. 


pressure  in  the  oxidant  gas  was  considered  to  be  increased  from  the 
inlet  to  outlet  of  oxidant,  thus  the  amount  of  transferred  water  from 
the  cathode  GDL  to  the  oxidant  gas  was  considered  to  be  decreased 
through  the  serpentine  channel  for  the  outlet.  The  evaporation 
margin  has  plus  value  in  the  entire  area  of  the  cell  in  Fig.  5,  and  it 
was  not  considered  that  the  flooding  occurred. 

The  current  density  distribution,  cell  temperature  distribution, 
water  evaporation  rate  distribution  on  cathode  GDL,  and  water 
evaporation  margin  distribution  on  cathode  GDL  for  the  conditions 
of  1  atm  and  counter-flow  are  shown  in  Figs.  6-9,  respectively. 


The  current  density  was  the  largest  at  the  oxidant  inlet  and  fuel 
outlet,  and  it  decreased  16%  through  the  serpentine  channel  for  the 
oxidant  outlet  and  fuel  inlet  in  Fig.  6.  The  cell  temperature  was  the 
lowest  at  the  oxidant  inlet  and  fuel  outlet,  and  it  increased  0.70% 
through  the  serpentine  channel  for  the  oxidant  outlet  and  fuel  inlet 
in  Fig.  7.  The  apparent  water  evaporation  rate  was  the  largest  at  the 
oxidant  inlet  and  fuel  outlet,  and  it  decreased  18%  through  the 
serpentine  channel  for  the  oxidant  outlet  and  fuel  inlet  in  Fig.  8.  The 
evaporation  margin  has  plus  value  in  the  entire  area  of  cathode  GDL 
in  Fig.  9,  and  it  was  not  considered  that  the  flooding  occurred. 


Fuel 
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Fig.  4.  Apparent  H20  evaporation  rate  (mol  s  1  cm  2),  1  atm,  co-flow. 
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Fig.  6.  Current  density  (A  cm  2),  1  atm,  counter-flow. 
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Fig.  7.  Cell  temperature  (K),  1  atm,  counter-flow. 
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Fig.  9.  H20  evaporation  margin  (mol  s  1  cm  2),  1  atm,  counter-flow. 


The  current  density  distribution,  cell  temperature  distribution, 
water  evaporation  rate  distribution  on  cathode  GDL,  and  water 
evaporation  margin  distribution  on  cathode  GDL  for  the  conditions 
of  3  atm  and  co-flow  are  shown  in  Figs.  10-13,  respectively. 

The  current  density  was  the  largest  at  the  inlet  of  oxidant  and 
fuel,  and  it  decreased  16%  through  the  serpentine  channel  for  the 
outlet  of  oxidant  and  fuel  in  Fig.  10.  Flere,  the  current  density 
distribution  was  similar  with  that  for  the  pressure  condition  of 
1  atm.  The  cell  temperature  was  the  lowest  at  the  inlet  of  oxidant 
and  fuel,  and  it  increased  0.99%  through  the  serpentine  channel  for 


the  outlet  of  oxidant  and  fuel  in  Fig.  11.  Flere,  the  cell  temperature 
distribution  was  lower  about  1  K  comparing  with  that  for  the 
pressure  condition  of  1  atm.  The  apparent  water  evaporation  rate 
was  lesser  significantly  at  the  outlet  of  oxidant  and  fuel  comparing 
with  that  for  1  atm  pressure  condition,  because  the  apparent  water 
evaporation  rate  was  decreased  according  to  the  increase  of  pres¬ 
sure  as  shown  in  equation  (2).  The  evaporation  margin  was  plus 
value  in  the  most  cathode  GDL  area  except  the  outlet  of  oxidant  and 
fuel.  At  the  outlet  of  oxidant  and  fuel,  the  evaporation  margin  was 
minus  value  and  it  was  considered  that  the  flooding  occurred. 
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Fig.  11.  Cell  temperature  (K),  3  atm,  co-flow. 
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Fig.  13.  H20  evaporation  margin  (mol  s  1  cm  2),  3  atm,  co-flow. 


The  flooding  area  corresponded  with  the  area  where  the 
apparent  water  evaporation  rate  was  less,  the  cell  temperature  was 
high,  and  the  current  density  was  less.  Here,  it  was  understood  that 
the  flooding  was  caused  by  the  low  apparent  water  evaporation 
rate.  In  addition,  the  flooding  was  considered  to  be  effected  by  the 
water  evaporation  and  condensation  with  oxidant  gas  flow, 
because  the  flooding  was  occurred  at  the  area  where  the  current 
density  was  low  and  the  amount  of  migrated  water  from  the  anode 
to  cathode  was  less. 


The  current  density  distribution,  cell  temperature  distribution, 
water  evaporation  rate  distribution  on  cathode  GDL,  and  water 
evaporation  margin  distribution  on  cathode  GDL  for  the  condition 
of  3  atm  and  counter-flow  are  shown  in  Figs.  14-17,  respectively. 

The  current  density  was  the  largest  at  the  oxidant  inlet  and  fuel 
outlet,  and  decreased  16%  through  the  serpentine  channel  for  the 
oxidant  outlet  and  fuel  inlet,  and  the  current  density  distribution 
was  similar  with  that  for  the  case  of  1  atm.  The  cell  temperature 
was  the  lowest  at  the  oxidant  inlet  and  fuel  outlet,  and  increased 
0.99%  through  the  serpentine  channel  for  the  oxidant  outlet  and 
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Fuel 


Fig.  15.  Cell  temperature  (K),  3  atm,  counter-flow. 


Fuel 


Fig.  17.  H20  evaporation  margin  (mol  s  1  cm  2),  3  atm,  counter-flow. 


fuel  inlet.  Here,  the  cell  temperature  distribution  was  lower  about 
1  K  comparing  with  that  for  the  pressure  condition  of  1  atm.  The 
apparent  water  evaporation  rate  was  lesser  significantly  at  the 
oxidant  outlet  and  fuel  inlet  comparing  with  that  for  1  atm  pressure 
condition.  The  evaporation  margin  was  plus  value  in  the  most  area 
of  cathode  GDL  except  the  oxidant  outlet  and  fuel  inlet.  At  the 
oxidant  outlet  and  fuel  inlet,  the  evaporation  margin  was  minus 
value  and  it  was  considered  that  the  flooding  occurred. 

As  similar  results  for  the  case  of  co-flow,  the  flooding  area  cor¬ 
responded  with  the  area  where  the  apparent  water  evaporation 


Fuel 


Fuel 


Fig.  16.  Apparent  H20  evaporation  rate  (mol  s  1  cm  2),  3  atm,  counter-flow. 


rate  was  less,  the  cell  temperature  was  high,  and  the  current 
density  was  less. 

The  flooding  areas  for  the  co-flow  and  counter-flow  of  3  atm 
pressure  condition  were  limited  at  the  oxidant  outlet  small  areas, 
and  the  areas  corresponded  to  the  areas  where  the  current  densi¬ 
ties  were  least.  Therefore,  the  effect  of  flooding  on  the  cell  perfor¬ 
mance  was  not  considered  to  be  significant. 

5.  Conclusion 

The  water  migration  from  anode  to  cathode  by  electrode  reac¬ 
tion,  the  water  transfer  by  evaporation  from  cathode  to  oxidant  gas 
with  oxidant  gas  flow  were  calculated  for  a  direct  humidify  model 
cell.  The  co-flow  and  counter-flow  for  the  fuel  and  oxidant  were 
considered  for  the  pressure  conditions  of  1  atm  and  3  atm,  and  as 
a  result  of  those,  the  cathode  flooding  characteristics  for  each 
condition  were  evaluated. 

In  the  case  of  pressure  condition  1  atm,  the  water  discharge  by 
oxidant  gas  was  conducted  sufficiently  and  the  flooding  was  not 
occurred  for  the  co-flow  and  counter-flow. 

In  the  case  of  pressure  condition  3  atm,  the  water  discharge  by 
oxidant  gas  was  prevented  at  the  outlet  of  oxidant  and  the 
flooding  was  occurred  for  the  co-flow  and  counter-flow  both  at 
the  oxidant  outlet.  The  flooding  area  corresponded  with  the  area 
where  the  apparent  water  evaporation  rate  was  less,  cell 
temperature  was  high,  and  current  density  was  less.  In  addition, 
the  flooding  was  considered  to  be  effected  by  the  water  evapo¬ 
ration  and  condensation  with  oxidant  gas  flow,  because  the 
flooding  was  occurred  at  the  area  where  the  current  density  was 
low  and  the  amount  of  migrated  water  from  the  anode  to  cathode 
was  less. 

The  flooding  areas  for  the  co-flow  and  counter-flow  of  3  atm 
pressure  condition  were  limited  at  the  oxidant  outlet  small  areas, 
and  the  areas  corresponded  to  the  areas  where  the  current  densi¬ 
ties  were  least.  Therefore,  the  effect  of  flooding  on  the  cell  perfor¬ 
mance  was  not  considered  to  be  significant. 
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List  of  symbols  O2  oxygen 

H20  water 


V  cell  voltage  (V) 

T  absolute  temperature  (I<) 

x  molar  fraction 

P  pressure  (atm) 

s  cell  area  (cm2) 

i  current  density  (A  cm'2) 

c  molar  concentration 

Ph2o  saturated  vapor  pressure  of  water  (atm) 

Qo  oxidant  molar  flow  rate  (mol  s-1) 

A  area  (cm2) 

Cm  Evaporation  margin  (mol  s_1  cm-2) 

Cc^o  H20  transfer  from  cathode  to  oxidant  (mol  s-1  cm-2) 
Ca^c  H20  transfer  from  anode  to  cathode  (mol  s_1  cm-2) 

Subscripts 
H2  hydrogen 
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